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Abstract Identifying the processes that determine the initial mass function of stars (IMF) 
is a fundamental problem in star formation theory. One of the major uncertain- 
ties is the exact chemical state of the star forming gas and its influence on the 
dynamical evolution. Most simulations of star forming clusters use an isother- 
mal equation of state (EOS). We address these issues and study the effect of a 
piecewise polytropic EOS on the formation of stellar clusters in turbulent, self- 
gravitating molecular clouds using three-dimensional, smoothed particle hydro- 
dynamics simulations. In these simulations stars form via a process we call 
gravoturbulent fragmentation, i.e., gravitational fragmentation of turbulent gas. 
To approximate the results of published predictions of the thermal behavior of 
collapsing clouds, we increase the polytropic exponent y from 0.7 to 1. 1 at some 
chosen density « c , which we vary from from 4.3 x 10 4 cirr 3 to 4.3 x 10 7 cm 4 . 
The change of thermodynamic state at n c selects a characteristic mass scale for 
fragmentation M c h, which we relate to the peak of the observed IMF. We find 
a relation M ch oc n~ 05±01 . Our investigation supports the idea that the dis- 
tribution of stellar masses largely depends on the thermodynamic state of the 
star-forming gas. The thermodynamic state of interstellar gas is a result of the 
balance between heating and cooling processes, which in turn are determined by 
fundamental atomic and molecular physics and by chemical abundances. Given 
the abundances, the derivation of a characteristic stellar mass may thus be based 
on universal quantities and constants. 
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1. Introduction 

Although the IMF has been derived from vastly different regions, from the 
solar vicinity to dense clusters of newly forme d stars, the ba sic features seem to 
be strikingly universal to all determinations (Kro upal 2001 ). Initial conditions 



in star forming regions can vary considerably. If the IMF depends on the initial 
conditions, there would thus be no reason for it to be universal. Therefore a 
derivation of the characteristic stellar mass that is based on fundamental atomic 
and molecular physics would be more consistent. 

There are many ways to approach the formation of stars and star clusters 
from a theoretical point of view. In particular models that connect stellar 
birth to the turbulent motions ubiquiteously observed in Galactic molecular 
clouds ha ve become in creasing ly popular in recent year s. See, e.g., the re- 
views bv Larsonl (J2003) and EVIac Low & Klessen (2004). The interplay be- 
tween turbulent motion and self-gravity of the cloud leads to a process we call 
gravoturbulent fragmentation: Supersonic turbulence generates strong density 
fluct uations with gravity taking over in the densest and most massive regions 
fe.g.lLarsorll 98ll:lFleckll 982l:IPadoanll 99.4IPadoan et alii 997l : lKlessen et all 
11 998L l2000l: lKlessenll200ll : IPadoan & Nordlundll2002h . Once gas clumps be- 
come gravitationally unstable, collapse sets in. The central density increases 
and soon individual or whole clusters of protostellar objects form and grow in 
mass via accretion from their infalling envelopes. 

However, most current results are based on models that do not treat thermal 
physics in detail. Typically, a simple isothermal equation of state (EOS) is 
used. The true nature of the EOS, thus, remains a major theoretical problem in 
understandi ng the fragment ation properties of molecular clouds. 

Recently iLietalJ (2003) conducted a systematic study of the effects of a 
varying polytropic exponent y on gravoturbulent fragmentation. Their results 
showed that y determines how strongly self-gravitating gas fragments. They 
found that the degree of fragmentation decreases with increasing polytropic 
exponent y in the range 0.2 < y < 1.4 although the total amount of mass 
in collapsed cores appears to remain roughly consistent through this range. 
These findings suggest that the IMF might be quite sensitive to the thermal 
physics. However in their computations, y was left strictly constant in each 
case. Here we study the effects of using a piecewise polytropic equation of 
state and investigate if a change in y determines the characteristic mass of the 
gas clump spectrum and thus, possibly, the turn-over mass of the IMF. 

2. Thermal Properties of Star-Forming Clouds 

Observational evidence predicts that dense p restellar cloud cores show rough 
balan ce between gravity and thermal pressure (Be nson & Mversll989l : lMvers et all 



1991). Thus, the thermodynamical properties of the gas play an important role 
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in determining how dense star-forming regions in molecular clouds collapse 
and fragment. Observational and theoretical studies of the thermal properties 
of collapsing clouds both indicate that at densities below 10~ 18 gem -3 , roughly 
corresponding to a number density of n = 2.5 x 10 5 cm~ 3 , the temperature 
decreases with increasing density. This is due to the strong dependence of 
molecular cooling rates on density ( Kovama & InutsukalEo OO). Therefore, the 



polytropic exponent y is below unity in this density regime. At densities above 
10" 18 g cm -3 , the gas becomes thermally coupled to the dust grains, which then 
control the temperature by far-infrared thermal emission. The balance between 
compressional heating and thermal cooling by dust causes the temperature to 
increase again slowly with increasing density. Thus the temperature-density 
relation can be a pproximated with y above unity i n this regime (Larson this 
volume; see also lLarsonlll985l : ISpaans & SilkEb oO). Changing y from a value 
below unity to a value ab ove unity results in a minimum temperature at the crit- 
ical density. As shown bv lLi et all (12003). gas fragments efficiently for y < 1.0 
and less efficiently for higher y. Thus, the Jeans mass at the critical density 
defines a characteristic mass for fragmentation, which may be related to the 
peak of the IMF. 

3. Numerical Approach 

To gain insight into how molecular cloud fragmentation the characteris- 
tic stellar mass may depend on the critical density we perform a series of 
smoothed particle hydrodynamics calculations of the gravitational fragmenta- 
tion of supersonically turbu lent molecul ar clouds using the parallel code GAD- 



GET designed by Springe! et all d200ll) . SPH is a Lagrangian method, where 



the fluid is represented by an ensemble of particles, and flow quantiti es are 
obtain ed b y averaging over a n appropriate subset of SPH particles, see Bend 



dl990h and lMonagharj dl992h . The method is able to resolve large density con- 
trasts as particles are free to move, and so na turally the partic l e con centration 
increases in high-density regions. We use the Bate & Burkert ( 19971) criterion 
the resolution limit of our calculations. It is adequate for the problem consid- 
ered here, where we follow the evolution of highly nonlinear density fluctu- 
ations created by supersonic turbulence. We repl ace the centra l high-density 
regions of collapsing gas cores by sink particles (IB ate et all 19951) . These parti- 
cles have the ability to accrete gas from their surroundings while keeping track 
of mass and momentum. This enables us to follow the dynamical evolution of 
the system over many local free-fall timescales. 

We compute models where the polytropic exponent changes from y = 0.7 to 
y - 1.1 for critical densities in the range 4.3 x 10 4 cm -3 <n c < 4.3 x 10 7 cm -3 . 
Each simulation starts with a uniform density distribution, and turbulence is 
driven on large scales, with wave numbers k in the range 1 < k < 2. We 
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Figure 1. Mass spectra of protostellar cores for four models with critical densities in the 
range 4.3 x 10 4 cnT 3 < n c < 4.3 x 10 7 cm -3 . We show two phases of evolution, when about 
10% and 30% of the mass has been accreted onto protostars. The vertical soli d line shows the 
median mass of the distribution. The dotted line serves as a reference to the S alpeteJ ll955l) 
slope of the observed IMF at high masses. The dashed line indicates our mass resolution limit. 



use the same driving field in all fo ur models. The global free-fall timescale is 
Tff w 10 5 yr. For further details see Uappsen et al.l (82004). 
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Figure 2. Plot of the median mass of the protostellar cores M me ^ m versus critical density n c . 
We display results for different ratios of accreted gas mass to total gas mass Af acc /M t01 , and fit 
the data with straight lines. Their slopes take the values -0.43 ± 0.05 (solid line), -0.52 ± 0.06 
(dashed-dotted line), and -0.60 ± 0.07 {dashed line), respectively. 



4. Dependency of the Characteristic Mass 

To illustrate the effects of varying the critical density, we plot in Fig. ^ the 
resulting mass spectra at different times when the fraction of of mass accumu- 
lated in protostellar objects has reached approximately 10% and 30%. This 
range of efficiencies co rresponds roughly t o the one expected for regions of 
clustered star formation (lLada & Ladal2003h . In the top-row model, the change 
in y occurs below the initial mean density It shows a flat distribution with only 
few, but massive cores. These reach masses up to 10 M Q and the minimum 
mass is about 0.3 M Q . The mass spectrum becomes more peaked for higher n c 
and shifts to lower masses. 

We find closest co rrespondence with the observed IMF ( Scaloll998 ; Kroupal 
120021: IChabrierll2003h for a critical density n c of 4.3 x 10 6 cm -3 and for stages 
of acc retion around 30%. For high masses, the distribution exhibits a (Salpeter 
power-law behavior. For masses about the median mass the distribution 
has a small plateau and then falls off towards smaller masses. 

The change of median mass M me( a an with critical density n c is quantified in 
Fig. El As n c increases M me di a n decreases. We fit our data with straight lines. 
The slopes take values between -0.4 and -0.6. 
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5. Discussion and Summary 

Using SPH simulations we investigate the influence of a piecewise poly- 
tropic EOS on the gravoturbulent fragmentation of molecular clouds. We study 
the case where the polytropic index y changes from 0.7 to 1.1 at a critical den- 
sity n c , and consider the range 4.3 x 10 4 cm -3 < n c < 4.3 x 10 7 cm -3 . 

A simple scaling argumen t based on the Jean s mass Mj at the critical density 
n c leads to Mj oc «~ u y5 (see Jappsen et al. If there is a close relation be- 



tween the average Jeans mass and the gravoturbulent fragmentation spectrum, 
a similar relation should hold for the characteristic mass M c h of protostellar 
cores. Our simulations qualitatively support this hypothesis, however, with the 
weaker density dependency M c ^ oc «~ a5±01 . So indeed, the density at which y 
changes from below unity to above unity defines a preferred mass scale. Con- 
sequently, the peak of the resulting mass spectrum decreases with increasing 
critical density. The distribution not only shows a pronounced maximum but 
also a power-law tail towards higher masses, similar to the observed IMF. 

Altogether, supersonic turbulence in self-gravitating molecular gas gener- 
ates a complex network of interacting filaments. The overall density distribu- 
tion is highly inhomogeneous. Turbulent compression sweeps up gas in some 
parts of the cloud, but other regions become rarefied. The fragmentation be- 
havior of the cloud and its ability to form stars depend on the EOS. However, 
once collapse sets in, the final mass of a fragment depends not only on the 
local Jeans criterion, but also on additional processes. For example, protostars 
grow in mass by accretion from their surrounding material. In turbulent clouds 
the properties of the gas reservoir are continuously changing. And in addition 
protostars may interact with each other, leading to ejection or mass exchange. 
These dynamical factors modify the resulting mass spectrum, and may explain 
why the characteristic stellar mass depends on the EOS more weakly than ex- 
pected from simple Jeans-mass scaling arguments. 

Our investigation supports the idea that the distribution of stellar masses 
depends, at least in part, on the thermodynamic state of the star-forming gas. 
If there is a low-density regime in molecular clouds where the temperature 
T sinks with increasing density p, followed by a higher-density phase where 
T increases with p, fragmentation seems likely to be favored at the transition 
density where the temperature reaches a minimum. This defines a character- 
istic mass scale. The thermodynamic state of interstellar gas is a result of the 
balance between heating and cooling processes, which in turn are determined 
by fundamental atomic and molecular physics and by chemical abundances. 
The theoretical derivation of a characteristic stellar mass may thus be based on 
quantities and constants that depend mostly on the chemical abundances in the 
star-lpwh klforming cloud. 
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